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We propose a scheme to dynamically realize a quantum memory based on the toric code. The code is generated 
from qubit systems with typical two-body interactions (Ising, XY, Heisenberg) using periodic, NMR-hke, pulse 
sequences. It allows one to encode the logical qubits without measurements and to protect them dynamically 
against the time evolution of the physical qubits. A weakly coupled cavity mode mediates a long-range attractive 
interaction between the stabilizer operators of the toric code, thereby suppressing the creation of thermal anyons. 
This significantly increases the lifetime of the memory compared to the code with nonintcracting stabilizers. 
We investigate how the fidelity, with which the toric code is realized, depends on the period length T of the 
pulse sequence and the magnitude of possible pulse errors. We derive an optimal period Tppt that maximizes the 
fidehty. 

PACS numbers: 03.67.Lx, 03.67 .Pp, 03.67.Ac 



I. INTRODUCTION 

One of the most promising proposals for the realization of 
a thermally stable quantum memory is based on topologically 
ordered phases of matter like Kitaev's toric code [1-3]. 

While the toric code allows for topological protection 
against local imperfections at zero temperature, it is suscep- 
tible to thermal fluctuations [4-7]: Anyons that are created at 
constant energy cost and move freely across the memory with- 
out additional energy penalty destroy the stored quantum in- 
formation in a time that does not increase with the linear size 
L of the memory. In fact, all two-dimensional (2D) and broad 
classes of three-dimensional (3D) stabilizer Hamiltonians with 
local interactions are subject to no-go theorems forbidding sta- 
biUty at finite temperature [8-10]. 3D stabilizer models with 
local interactions that do not satisfy the criteria of applicability 
of the no-go theorems have been proposed and studied. None 
of them, however, is so far expected to be stable at any nonzero 
temperature [11-14]. 

Despite the intrinsic thermal fragility of the toric code, it is 
possible to considerably improve the memory lifetime by al- 
lowing for long-range repulsive interactions between anyons 
as proposed in Ref. [15] and later studied in several systems in 
Refs. [16-18]. In fact, this leads to a suppression of the anyon 
density in the thermodynamic limit and thus to a memory 
lifetime increasing polynomially with L. Such memories are 
called self-correcting, as their stability against errors caused 
by the thermal environment is "built in" in the sense that no 
active (measurement-based) error correction is required. In 
similar approaches, long-range attractive interaction between 
anyons as a way to suppress their motion across the memory 
has been proposed in Ref. [2] and later studied in Ref. [19] by 
coupling the toric code to a bosonic bath. In this case the diffu- 
sion of anyons is reduced by the attractive interaction and the 
memory lifetime increases polynomially with L, but the model 
requires unbounded-strength interactions between anyon oper- 
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ators and the bosonic bath [19]. Recently, a three-dimensional 
model where toric code stabilizers are locally coupled to the 
spins of a ferromagnet has been proposed in Ref. [20]. The 
attractive interaction between stabilizers is then mediated by 
Goldstone modes and leads to a memory lifetime increasing 
exponentially with L. 

In this work, we study a similar model, for which the at- 
tractive interaction between stabilizers of the toric code is me- 
diated by cavity modes. This model was first proposed in 
Ref. [15] and studied in detail in Ref. [16] as a low-energy ef- 
fective Hamiltonian of an anisotropic honeycomb model cou- 
pled to a cavity mode. In the same sense as for these cavity- 
based setups, the model proposed here implements a passive 
error correction which renders it self-correcting against ther- 
mal fluctuations. This self-correcting property, however, is 
limited by pulse imperfections as weU as the validity of the 
perturbation theory used to derive the efl'ective stabilizer in- 
teraction. 

The toric code Hamiltonian contains many-body interac- 
tions that are not directly realized in nature. Nevertheless, sev- 
eral ways have been proposed as to how the toric code (and 
similar stabilizer codes) can be implemented in practice. Be- 
sides the low-energy limit of the honeycomb lattice [21, 22], 
it can emerge dynamically as a coarse-grained description of 
a quantum simulation with discrete time steps using Rydberg 
atoms [23], ions [24, 25], or polar molecules [26, 27] in op- 
tical traps. This work is based on another kind of dynamical 
implementation that employs NMR-like, periodic sequences 
of short external pulses to induce the dynamics of the code 
Hamiltonian in solid-state systems, similar to a recent proposal 
in Ref. [28]. By its use of periodic pulses, this scheme of dy- 
namically generating a desired Hamiltonian dynamics is re- 
lated to the so-called dynamical decoupling methods [29, 30] 
that mainly aim to suppress the effect of a given or even un- 
known system-bath interaction on (time-)average. From an 
abstract, mathematical point of view, both kinds of pulsing 
schemes are rather similar insofar as their effect can be cap- 
tured by taking the time average of Hamiltonian terms or error 
generators that are transformed by appropriate unitary oper- 
ations (the pulses). Compared to the dynamical decoupling 



schemes, the method used here is particularly suited to gener- 
ate the Hamiltonian dynamics of the planar code with as few 
linewise applied rotations as possible. Also, in addition to the 
suppression of decoherence that is caused by a decouplinglike 
effect of the pulses, our proposal achieves an even stronger sta- 
bilization against the thermal environment by a coupling to a 
nonlocal field (the cavity mode). 

Starting from a system of 1/2 spins coupled to a cavity 
mode, we show how to generate a toric code with long-range 
attractive interactions between stabilizers. For perfect pulses, 
the fidelity F with which the dynamics of the code Hamilto- 
nian is generated only depends on the structure and time period 
T of the generating sequence reaching the theoretical limit of 
F = 1 for T — > 0. In reality, however, unavoidable pulse im- 
perfections require to minimize the number of pulses per time 
and the fidelity reaches its maximum for a finite period Topt. 
We derive an explicit formula for Topt and show that for pulse 
errors of a relative magnitude around 10^* a fidelity of more 
than 0.99 can be realized over a long time t > lOOA^^, where 
A is the energy scale to create a bare anyon. Finally, we show 
how to encode a logical qubit by a sequence of magnetic pulses 
and thus without the need of stabilizer measurements. 

The paper is organized as follows. After the model system 
is introduced in Sec. II, we show how to generate the dynamics 
of the toric code Hamiltonian with a periodic pulse sequence 
in Sec. III. We study both numerically and analytically how 
the gate fidelity F of the induced time evolution depends (i) 
on the structure of the pulse sequence and (ii) on the magni- 
tude of random pulse errors in Sec. IV. Furthermore, based 
on both these limiting factors for F, an analytic expression 
for the optimal sequence period length Topt > is derived 
that maximizes the fidelity. In Sec. V, we derive an effective 
low-energy Hamiltonian for a pulsed system of qubits that are 
weakly coupled to a cavity mode. The cavity induces an at- 
tractive stabilizer interaction that protects the system against 
errors caused by a thermal environment. We explain that, with 
respect to these errors, the resulting memory can be consid- 
ered as self-correcting. Finally, the measurement-free encod- 
ing of logical qubit states into the ground state manifold is ex- 
plained in Sec. VI. Appendices A-C contain details about the 
pulse operations for arbitrarily large qubit arrays and technical 
derivations regarding the fidelity dependence on the sequence 
structure and on the magnitude of pulse errors. 

n. THE MODEL SYSTEM 

We consider a quadratic lattice of noninteracting qubits 
with site-independent level splitting, which are weakly cou- 
pled (0 < (5 <C A) to a single cavity mode 

HoiS) = -4[A + ,5(6 + b^)] ^ ai^') + huob^b, (1) 

where cr^'' with >c G {x, y, z} denotes a Pauli matrix acting 
on the qubit on lattice site j. Operator b annihilates a photon 
of energy huo in the cavity. We assume that two-qubit gates 
can be applied by an external switching on and off of a "nat- 
ural" (system-dependent) two-qubit interaction. For example. 




FIG. 1. (Color online) Schematic of the planar code Hamiltonian 
equation (2) on the square lattice (dotted lines). Plaquettes (P) are 
associated with the edges of a unit cell and stars (S) with the edges 
of the cells around a vertex. A circle indicates a Pauli operator act- 
ing on a qubit and the wiggly lines connect operators that have to be 
multiplied to obtain the corresponding stabilizer operator. 

in the case of spin qubits realized in single-electron quantum 
dots an effective nearest-neighbor Heisenberg interaction can 
be switched electrically by changing the transparency of the 
tunneling barrier between two dots using gate electrodes [31]. 
The goal of this work is to show that with the choice of a 
proper periodic sequence of single-qubit rotations and two- 
qubit gates, a self-correcting topological quantum memory can 
be generated dynamically. We show this for the case of the pla- 
nar code (toric code with boundaries), which is described by 
Hamiltonian 

Hpc^-AfY^Ap+Y."^^) --^Y.^- ^2) 

^ p s ^ a 

For a given quadratic lattice of qubits sitting on the edges of a 
unit cell, indices p and s run over all plaquettes and stars, re- 
spectively. A plaquette is the set {pi, ... ,^4} of qubits sitting 
on the edge of a single unit cell and the corresponding operator 
is given by Ap = Y[i o'i^'''. Associated to every star is the set 
{si, . . . , S4} of qubits around a vertex connecting four neigh- 
boring cells, where As — Yli '^i''^^ ■ At the edges of the code, 
the Ap and As consist of only three Pauli operators. Figure 1 
illustrates the structure of Hpc schematically. The stabilizer 
operators Wa are identical to the plaquettes and stars and in- 
troduced to obtain a simpler notation. Hence, index a runs 
through all unit cells and vertices of the quadratic lattice. 

We want to stress, however, that the method described here 
is not restricted to the planar code or the particular free Hamil- 
tonian Hq. In fact, it can in principle be used to generate a 
quite large class of code Hamiltonians. Restrictions of this 
method rather concern the practical implementation of a par- 
ticular scheme in real experiments. Details such as the re- 
quired fidelities on certain time scales, addressability issues, 
the accuracy of external operations, etc., have to be examined 
with respect to a particular physical system. Below we present 
an example of how to implement the dynamically generated 
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topologically protected quantum memory based on Hpc in re- 
alistic systems. 



ni. DYNAMIC GENERATION OF PLANAR CODE 
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In this section, we explain in detail how external pulses and 
two-qubit operations can be used to dynamically "generate" 
the planar code Hamiltonian from the qubit part of ffo = 0) . 
Here and in the following "generating a Hamiltonian -ffav" 
is used to mean that, due to external pulses applied between 
certain discrete (stroboscopic) times, a system with Hamilto- 
nian H evolves as if its dynamics was governed by i?av — the 
"average Hamiltonian." After a short introduction to average 
Hamiltonian theory [32], we proceed by showing one possible 
sequence of pulses that generates H^c from Ha{5 — 0). Fi- 
nally, a nonzero coupling of the qubits to a cavity mode [as in 
Eq. (1)] can be utilized to induce a strongly nonlocal interac- 
tion between the anyon operators of H^c- 
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FIG. 2. (a) Generic pulse (operation) sequence. The unitary oper- 
ations Ri (gray columns) with YYi=i Ri = interrupt periods of 
(free) propagation with Hq. (b) Structure of a sequence of duration 
T to generate a plaquette or star operator. A sequence of the same 
structure is employed to generate (each) one of the quarters of the 
planar code. 



A. Average Hamiltonian Tlieory 

Average Hamiltonian theory describes how time-periodic 
(externally controlled) unitary transformations can be used to 
let the evolution of one system mimic that of another system of 
the experimenter's choosing. For the purpose of our paper, it is 
sufficient to consider the case of n-|- 1 periods to, ... ,tn of free 
propagation with Hq, which are separated in time by n unitary 
transformations Ri with i = 1, .... r? and Yli=i Ri = 
is assumed, that all the Ri can be performed within a typical 
time r, , where the Ti define the smallest time scale of the sys- 
tem. This is illustrated in Fig. 2(a). The physical time of the 
entire sequence is then given by T = to + J2i=ii^i + '''«)' 
after which the model will have evolved according to the time 
evolution operator Ut = Uo(tn)Rn ■ ■ ■ R2Uoiti)RiUo{to) = 
exp{— jfifav} with f = X;r=o *^ ^ T. Here and in the fol- 
lowing, h is set to 1 and we defined Uo{t) = cxp{— it_ffo}- 
With the Magnus expansion [32, 33], such a product of uni- 
tary operators can always be written as a single exponential 
of the average Hamiltonian ffav Hence, if the pulse sequence 
is applied periodically, at integermultiples ofT (stroboscopic 
times) the model system will have evolved as if governed by 
Hamiltonian THay/T. In general, the exact ffav can only be 
given as an infinite expansion in orders of TA, where is 
(of the order of) the characteristic time scale of the unpulsed 
system Hq: 

ffav = ffi°^+ffW+ffi,2^.. (3) 

For small enough TA, however, it is often sufficient to orily 
consider the lowest-order term 

1 " 

H..^Hi^^ = -Y,H,ti, (4) 
with Hj = 7^]ffo7^J for j ^ and % = W^^^ Rk. 



B. Generating Sequence for Planar Code 

In order for a sequence as shown in Fig. 2(a) to be usable 
in real experiments and applications it has to be as short as 
possible and the operations Ri have to be decomposable into 
simple elementary operations. We will show how ffpc can be 

generated from Ho{5 = 0) using only (i) 7r/2 single-qubit ro- 
tations about the x and y axes and (ii) the two-qubit controlled 

phase gate [/pG = cxp{i7rcrP ai^^ /4} = (1 +i f7i^V?V\/2 
between qubits on neighboring sites (denoted by 1 and 2). 
Each of these operations only have to be applicable on rows 
(columns, diagonals) of qubits simultaneously. In particular, 
we will not require that physical qubits can be addressed indi- 
vidually, which might be advantageous for certain qubit archi- 
tectures. 

In general, besides single-qubit rotations the generation of 

ffpc only requires an operation that allows one to transform 

(%>) 

Pauli terms h = a)i in the average Hamiltonian according to 
— >■ cr^^V^^) for at least one axis k. Here, p = 1, 2 and 
p ^ p are the site indices of neighboring qubits. The spin axes 
X, x', and do not need to be different. 

Gate UpG is an example of such an operation, which can 
be seen by considering a sequence of two pulses with TZi = 
Ri = UpG, and to = ^2 = so that ti = T. According to Eq. 
(4), we obtain the average Hamiltonian to be h,^y = UpQhUpc- 
A straightforward calculation shows that this transformation 
has nontrivial effects only on terms that contain exactly one 
of either ax or ay. Concretely, the transformation converts 

a'i^ ^ a^^^af^ and a^^^ ^ -af^af\ 

In a system, in which the qubit interaction is Ising-like with 
-ffising = —JJ2{jk) (^z'^f^z'^ and J > 0, where the sum in- 
cludes all next neighbors {jk), the operation Upc can be per- 
formed by switching the interaction on for a time r = 7r/(4J). 
Note that similar operations can also be realized based on 
other nearest-neighbor interactions such as XY or Heisen- 
berg. In the XY case, for example, this role can be played 
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by exp{±i Jr(cri^'cri^^ + a'^y^a^y^)} for r = 7r/(4J), while 
using a (sightly modified) pulse sequence of similar complex- 
ity to the one we present for the Ising case. Analogous opera- 
tions can be found in the Heisenberg case. Thus, our scheme 
is not restricted to qubit systems with a particular kind of two- 
particle interaction. 

We proceed by showing how to generate a plaquette or star 
operator for a 2 x 2 quadratic lattice of free qubits. For sim- 
plicity, we consider the case A = 1/4 and (5 = 0, given 

by H{) = "'z"'^' ^i'^h j denoting the qubits as in Fig. 

3(a), and _ffav — Ag = Y\- cfx^ The pulse sequence con- 
sists of four (complex) operations with — = Q and 
2ti — 2^3 = ^2 = T/2. Its structure is illustrated in Fig. 
2 (b), where we have to specify operations i?i and i?2 so that 
the sequence yields the desired (lowest-order) average Hamil- 
tonian. 

A possible sequence of operations to generate a fourth-order 
spin Hamiltonian starting from i/o is illustrated in Fig. 3, while 
Table I shows the terms of the transformed Hamiltonian for 
each step. After rotating the qubits along one diagonal [upward 
in Fig. 3(b)] about the y axis by 7r/2, second-order terms (indi- 
cated by wiggly lines) are generated by applying Upc (hatched 
areas) to two parallel edges of the lattice [vertical in Fig. 3(c)]. 
Since applying the entangling gate to the remaining (horizon- 
tal) edges in the resulting configuration would yield third- in- 
stead of fourth-order terms, the qubits along one of these edges 
[the lower in Fig. 3(d)] have to be rotated, so that the z and 
y terms are interchanged. Subsequent application of C/pc to 
the remaining edges [horizontal in Fig. 3(e)] now generates 
[amongst others, see Eq. (6)] a fourth-order term. 

The last step shown in Fig. 3(f) is not actually necessary in 
case of the simple 2x2 lattice. Depending on whether a pla- 
quette or staright-hand-sideto be generated, one could either y 
rotate the rightmost or leftmost edge after step (e). In view of 
the application to a larger planar code like in Fig. 1, however, 
we interchange the lower edge qubits a second time [Fig. 3(f)]. 
A final rotation of the qubits along a properly chosen diagonal 
then yields the desired fourth-order operator (not shown in the 
figure). In the case of the planar code on a larger lattice, the 
additional step is required due to the presence of (third-order) 
boundary operators in H^c- It arranges their components in 
straight lines with those of the interior operators, so that a final 
rotation along properly chosen diagonals generates a Hamilto- 
nian with only z and x terms, respectively (cf. Appendix A). 

The candidate for Ri to generate a star operator from is 
therefore the product of all the operations of steps (b)-(f) plus 
a y rotation of qubits 1 and 4 with 

n, „ rr(2,3) ,r(l,3) ..(2,4) ..(3,4) ..(1,2) ..(3,4) ..(3,4) ..(1,4) 

(5) 

where J/pG*^' entangles qubits j and k and Ux^z denotes a 7r/2 
rotation about the axis perpendicular to x and z. Note that 
Ri acts on the system by performing the elementary steps in 
reverse order. If we have only Ri or, equivalently, set R2 = 1., 
we obtain average Hamiltonian 

H, = + a^^ + af^a^^^ - a^^\ (6) 




FIG. 3. (Color online) First 6 of 7 elementary steps to generate fourth- 
order terms in an average Hamiltonian starting from a single-particle 
Hamiltonian for qubits on a quadratic lattice. 



TABLE I. Hamiltonian terms for each step in Fig. 3 (up 
to a sign). Pauli operators that changed compared to 
the previous step are shown against a gray background. 
A final y rotation (not shown) of qubits 1 and 4 yields 
Eq. (6).' 
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Apart from the desired star operator, it contains unwanted 
single- and two-particle by-products (see columns 1,3, and 4 in 
Table I). These can be removed from the average Hamiltonian 
with an R2 that yields a Hamiltonian H2 = Ri^ R^^ H0R2R1 
[see Eq. (4)], whose first- and second-order terms have the 
opposite sign compared to Hi. Note that if such an opera- 
tion i?2 were considered as just another step of the sequence 
in Fig. 3, it would have to be performed prior to (b): With 
T/q := R2^HqR2, we can write H2 = Ri^HqRi. Hence, we 
need to find an operation R2 that changes Hq to Hq in such a 
way that the subsequent application of Ri results in the desired 
sign flips. 

It turns out that the simple ansatz to just rotate qubits 1 
and 3 by an angle of tt about either the y or z axis (in an 
attempt to add one sign to each first- and second-order term, 
while adding two signs to the fourth-order operator) will not 
do the trick. A straightforward calculation shows that qubit 3 
has to be rotated as well to give the correct result. Thus, with 
i?2 ~ Uy^'^''^\Tr) and using — Hi we obtain the aver- 
age Hamiltonian H^y — {Hi + i?2)/2 — <ji^^ cr^x^^ a^'' 

Here, operator ui^^'"''^"\<p) rotates qubits ii,...,i„ about 
axis >c — x,y, z hy angle 4>. Finally, R2 has to be decom- 
posed into operations that are applied linewise. One possible 
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FIG. 4. Numerical gate fidelity versus time of a dynamically gener- 
ated planar code of = 13 qubits (inset) coupled to a cavity mode 
with three dimensional Fock space. For A = loq, 5/ A = 0.1, and 
T = A~^/8, already the lowest-order average Hamiltonian yields 
F > 99.9% for times t < lOOA"^ (solid line). Adding the second- 
order contribution decreases the deviation from perfect fidelity by an- 
other two orders of magnitude (dashed and dotted line, plot in inset). 
Based on the two leading orders of the average Hamiltonian, the ex- 
act fideUty can be approximated reasonably well, both numerically 
(dashed Une) and analytically (dotted Une). 

way to do this is 

i?2 = u^'''\n/2) c/i^'^) w u^'''K^m. (7) 

With a pulse sequence of the exact same structure and sim- 
ilar (generalized) i?i and R2, we can also generate either mul- 
tiple plaquette or star operators on an arbitrary large quadratic 
lattice. Details about both operations are given in Appendix 
A. We point out that only one-quarter (even or odd subsets 
of plaquettes or stars) of Hpc can be generated with one such 
sequence. This is due to the fact that only disjoint, i.e., half 
of the pairs of qubits, can be coupled in an entangling step per 
lattice dimension, to generate fourth- and third-order terms of 
the correct structure. Since each quarter of the code is gener- 
ated for one-quarter of the time only, the single-particle energy 
of Hamiltonian Hq (5 = 0) from Eq. (1) was set to be — 4A to 
yield a Hpc with energy gap — A. 

IV. FIDELITY LIMITS DUE TO SEQUENCE STRUCTURE 
AND PULSE ERRORS 

A suitable measure for the similarity of the dynamically 
generated dynamics with that of the toric code Hamiltonian 
is given by the gate fidelity 

F{t) = \Ti[exp{itH,,)Uj.{t)]\ /2^, (8) 

where Up{t) denotes the evolution operator of the pulsed sys- 
tem. The gate fidelity indicates how well an arbitrary basis 
of the whole Hilbert space evolves according to i?av By con- 
trast, a state-specific fidelity only quantifies how well the dy- 
namics is described by a particular H^y for that same specific 



set of states. For example, a high value for a fidelity that mea- 
sures how much of an encoded state remains in the ground 
state manifold over time, is consistent with the dynamics of 
any Hamiltonian that coincides with that of ifav on the ground 
state manifold. For the generation of the stabilizer interaction 
by means of a cavity, however, it is important that besides the 
ground states all excited states propagate according to H^y, as 
well. Hence, in contrast to state fidelities, only a gate fidelity 
F{t) » 1 for iJav = Hpc can be used as an indicator that the 
cavity indeed induces the stabilizer interaction. 

For decreasing sequence duration of the order of or shorter 
than A~^, by simply generating each four quarters alternat- 
ingly with a period of T/4, the lowest-order average Hamilto- 
nian 

Hi°^ = [1 + S/A{b + b^)]Hpc + ujob^b (9) 

describes the system dynamics with increasing fidelity over 
longer times even between the end points of the sequence. 
With a suitably symmetrized version (cf. Fig. 8 in Ap- 
pendix B), all odd-order contributions to H.^^ [see Eq. (3)] can 
be completely suppressed, so that the leading-order deviation 
is given by 

= "^QI^Q^ + + b^)] with 

3A -5 1 6 7 ^ 4 3' 
\-5 1 7 12/ \-5/ 

where Q = (Qi , • • • , Q4) is a vector consisting of the four 
quarters Qi = —A j T^a; of the planar code defined by 
appropriate, disjoint sets of anyon indices {oj}. 

Figure 4 shows the time-dependent gate fidelity for perfect 
pulses and a lattice of = 13 qubits (inset) coupled to a cav- 
ity with three-dimensional Fock space {|n ± 1), which 
is evaluated with a numerically exact Chebyshev expansion 
of the time evolution operator [34]. The system parameters 
are T = A'^/S, S = O.lA, and ojq = A. Akeady for 

ifav = ^^av^ (soUd Une), the fidelity does not drop below 
99.9% before t = lOOA"'^, while adding the second-order 
contribution from Eq. (10) decreases the deviation from the 
perfect fidelity by another two orders of magnitude (dashed 
and dotted line, plot in inset). In reality, however, the fidelity 
will be lower than this theoretical maximum due to errors in 
the pulsing scheme, decoherence by the noisy environment, 
and fluctuations of the microwave beam that is used to induce 
anyon interactions. Compared to systems that directly realize 
the toric code Hamiltonian (rather than just its dynamics), e.g., 
in some low-energy hmit of some suitable local lattice Hamil- 
tonian, these effects will decrease the Ufetime of codewords. 
Nevertheless, the numerical results indicate that "* + -ffiv ■* 
describes the perfect-pulse dynamics of the system very well 
and therefore allows one to estimate the size of deviations from 
the intended Hamiltonian (9) that arise solely by the structure 
of the pulse sequence. This is illustrated by the dashed line in 
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Fig. 4, which gives 

|Tr[exp(z7Ji«)i)exp(-z{i?i°) + 

(11) 

as resuh of a numerically exact calculation for the 1 3-qubit sys- 
tem. Its quadratic behavior is very similar to that of the pulsed 
system, while the quantitative deviations are due to terms of 
order 4 and higher In Appendix B it is shown that i^'^' can 
be approximated by F'--'^^ « 1 — Cav^^ with 

for planar codes of arbitrary linear size L, which we define as 
the length of the larger quadratic sublattice for codes arranged 
as in Figs. 1, 4, and 7. For example, the code with 13 = 3^ + 2^ 
qubits has length L ~ ?>, while the code with 25 = 4^ + 
3^ qubits has length L = 4. For L = 3, the approximate 
is given by the dotted line in Fig. 4. Although Eq. (12) 
systematically overestimates the fidelity of the pulsed system 
F, it is sufficiently accurate to provide the bounds 1 < ( 1 — 
F) /{cm't^ ) < 2 for deviations due to higher-order corrections 
to i?iv\ as long as c^yt^ <C 1. 

While the deviations that are due to the structure of the pulse 
sequence decrease rapidly for smaller periods T, errors caused 
by imperfect pulses will generally increase with the pulsing 
rate. As a consequence, a finite magnitude of pulse errors 
entails an optimal, finite value of T, for which the total de- 
viations due to both error sources are minimal. We estimate 
the effect of pulse imperfections using a simple model with 
errors manifesting in small random deviations (50 <C 1 from 
the intended angle Oq in R{d9) = cxp{i[dQ + S9]S/2), where 

(?) (1) (7+1) 

S — a)i for rotations about axis x and S = for 

the phase gate, respectively. The error distribution with width 
= \J {60"^) > is assumed to be Gaussian, equal for all 
considered kinds of pulses, and unbiased {{59) — 0). In Ap- 
pendix C, we show that this error model leads to an average 
gate fidelity (F) w 1 — Cen-t — ac^^t^ with 1 < a < 2 and 

2 

Ce„=(18L2-10L + 7)g; (13) 

for Ceni ^ 1- As shown in Fig. 5, this simple approximation 
(solid lines) agrees well with the numerical results (crosses) of 
the L = 3 code for T = A~ 78 and pulse errors < 10"^ x^e 
fidelity exhibits a linear and quadratic behavior in dependence 
on the time and on the magnitude of pulse errors, respectively. 
Thus, by minimizing the total deviations Cen-i + acav^^ with 
respect to T, we can estimate the optimal period length to be 

_ / 96[2L(9L-5) + 7]ag y/^ 
°'"~\[hL{L-l)5u:^l^fat) ■ ^ ' 

For example, to maximize (F) for the L = 3 code with 
ujQ = A, S ~ O.IA, and an error magnitude of ag = 0.001 
after a propagation time of t = lOOA^^, we have to set 
Topt « O.lA^^, which results in a fidelity of about 91%. For 
ag = 10-4, we obtain Top, w 0.04A-1 and (i^)(100A-i) > 
99%. When realizing the optimal sequence period, the total 



xlO" 




'^e t (units of T) 

FIG. 5. (Color online) Average of numerical gate fidelity (F) with 
imperfect pulses and T = A~^/8. (a) shows (F) — 1 versus the 
width ag of the Gaussian distributed pulse errors for t = A"^, while 
(b) shows the fidelity versus time for ere = 5 x 10^*. The solid lines 
are evaluated using Eq. (13) and agree well with the numerical data 
(crosses). 



deviation from perfect fidelity scales with the linear size of the 
memory like 0(Li2/5) This Hmits the gain in stability that is 
achievable by the induced stabilizer interaction (see Sec. V). 
In an experimental realization this leads to an optimal L that 
maximizes the lifetime of the memory. 



V. SUPPRESSION OF THERMAL ANYONS DUE TO 
CAVITY-INDUCED STABILIZER INTERACTIONS 

In this section, we derive the effective long-range attractive 
interaction between stabilizers that is induced by the coupling 
of a cavity mode to the entire lattice of physical qubits. As 
long as the coupling is weak (6 <^ ujq, A), a Schrieffer- Wolff 
transformation [35] yields the effective low-energy Hamilto- 

nian for a cavity that is kept in a Fock state \n) with a constant 

(2) 

number n of photons. Aside from small deviations due to H^^ 
and pulse imperfections, the total average Hamiltonian is then 
given by 

=-AVWa-— V WaW,'. (15) 

a a^a' 

The interaction term, which describes an effective long-range 
interaction between the stabilizer operators, couples every pair 
of stabilizers of the code via absorption and emission of a 
single photon. Note that this term in the effective Hamilto- 
nian (15) is at least two orders of magnitude stronger than the 
one that appears in the second-order average Hamiltonian in 
Eq. (10) for typical values of Top, and (5 ^ ~ A. 

It was shown in Ref. [16] that the energy penalty for the 
creation of anyons grows with for a code Hamiltonian (15) 
that features an attractive constant long-range interaction be- 
tween stabilizers. As a consequence, the lifetime of logical 
qubit states (or codewords) due to independent single-qubit er- 
rors that are created by the thermal environment increases ex- 
ponentially with or, equivalently, with the number of sta- 
bilizer operators Wa- In a system that is a direct realization 
of Hamiltonian (15), the density of anyons therefore goes to 
zero in the thermodynamic limit and without the need for any 
measurements-based, active error-correction procedure, the 
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memory retains its information indefinitely. In other words, 
the couphng to the cavity mode allows one to passively pro- 
tect the quantum information against thermal fluctuations. In 
this sense the memory is called self-correcting. 

Note that, while being an element of repeated external ma- 
nipulation, the pulsing of the system does not in itself consti- 
tute an active error-correction procedure: Neither does it re- 
quire one to extract information about the system state at any 
point in time (by measuring some system properties) in order 
to detect an error, nor to manipulate the system in a way that 
is conditional on the result of (such) an error syndrome mea- 
surement. Rather, any error-correcting effect is purely passive 
and achieved by employing a time invariant sequence of pulses 
to create an environment that hinders anyon creation. Conse- 
quently, we can regard the time invariant pulsing as an integral 
part of the system — a system that is passively protected against 
thermal fluctuations and therefore self-correcting in the sense 
explained above. In contrast to a direct realization of Hamil- 
tonian (15), however, the dynamical quantum memory has a 
finite lifetime in the thermodynamic limit. Besides the factors 
discussed in Sect. IV, it is limited by (i) fluctuations in the cav- 
ity mode, (ii) a breakdown of the perturbation theory used to 
derive the effective interaction in Eq. (15) for large L, and (iii) 
system-bath interaction processes with coherence times that 
are of the order of or shorter than T. 

In summary, the following constraints between energies and 
time scales have to be fulfilled by every implementation of the 
dynamically generated quantum memory: 

(5«A,c^o, r«Ti,r2, /3-i«A, (16) 

where f3~^ is the thermal energy and Ti and T2 are the relax- 
ation and dephasing times of the physical qubits, respectively. 
The latter provide time scales both for the rates of errors cre- 
ated in the code as well as for the process of error creation it- 
self Further restrictions may appear depending on the details 
of the physical realization. 

VI. PREPARATION OF CODEWORD STATES 

In addition to generating the dynamics of the planar code 
Hamiltonian, the pulsing method can also be used to prepare 
codewords {|0), |1)} without the need to perform measure- 
ments on the system [28]. Typically, a codeword of the pla- 
nar code is prepared by consecutive projective measurements 
of all (commuting) star operators performed on initial state 
|0, 0, . . . , 0) that is obtained by preparing all physical qubits 
in the spin-up state and which is already a simultaneous eigen- 
state of all the plaquette operators with eigenvalue(s) +1. 

Specifically, a codeword in the ground state multiplet of 
Hpc — one that is free of anyons — is obtained by 

s ^ 

The nonunitary projection operators c>c (1 + As) cannot, in 
general, be implemented with external pulses that are essen- 
tially unitary operations. However, if for one of the qubits 




step one step two 

FIG. 6. (Color online) (a) and (b): Two steps for preparing jO) from 
initial state |0, . . . , 0) for tiie L = 3 code. The modified star oper- 
ators are obtained by rotating one ax to —(Jy. In each of the steps, 
tiie indicated modified stabilizers are generated for a (effective) time 
t — 7r/(4A) and the ay are chosen such that none of the stabilizers 
from step one acted on y-qubits from step two. (c) Average codeword 
fidelity {F}c as a function of pulse errors. Numerics (crosses) agree 
well with the analytical result (solid line). 



Sfc G {si,---,S4}of star Ag, the state prior to the projec- 
tion can be written as = \0) ^ \(j>), where |0) is an 
arbitrary state in the joint Hilbert space of all — 1 remain- 
ing qubits, then, by dynamically generating —AAs{k) with 

As{k) = -a^""^ Hi^^fc '^^^'"^ for ^ time t = 7r/(4A), we get 
expM,(fc)/4]|7^,) = ^(1 + As)\A), (18) 

which can be traced back to the identity ~iay\0) — cra:|0). 

As per Eqs. (17) and (18), one As{k) for each star 
operatoright-hand-sideto be generated, to prepare state |0). 
The order of their application and the set of rotated qubit terms 
has to be chosen such that for each As{k), none of the previ- 
ously applied As' {k') has acted on the qubit s^. For the L ~ 3 
code, this is a two-step process and schematically illustrated in 
Fig. 6. In each step, one-half of the modified star operators are 
generated. While there are 4x3x3 equivalent possibilities to 
choose three spin operators to rotate from a; to ?/ in step one, the 
corresponding operators in step two are uniquely determined. 
This procedure works particularly well for the L — 3 code, as 
it mainly consists of edge operators; for L > 3, the preparation 
requires a larger number of steps. It is clear, however, that such 
a set of steps can always be found and generated with linewise 
rotations and entangling gates only: In the worst case, the pro- 
cedure consists of L{L — 1) steps — one per star operator. 

Figure 6(c) shows the codeword fidelity, i.e., the probability 
Fc ■■= |(0|C/prep|0, . . . ,0)1^ to find the system in State |0) after 
the preparation sequence J/prep, as a function of the pulse errors 
and with (5 = for the time of preparation [36]. Again the fi- 
delity decreases quadratically with ae and agrees well with the 
simple approximation from Appendix C: Counting the number 
of mutually inverse pulses in the preparation sequence corre- 
sponding to Figs. 6(a) and 6(b) yields Fq ~ 1 — 63(jg for the 
L = 3 code. 
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Vn. CONCLUSIONS 

We showed how to dynamically realize a quantum mem- 
ory based on noninteracting qubits that is stabilized against 
thermal fluctuations by a weakly coupled microwave cavity. 
Properly designed, periodic sequences of pulses, implement- 
ing single-qubit rotations and controlled phase gates, can be 
used to induce the time evolution of Kitaev's toric code Hamil- 
tonian. This allows one to prepjire logical qubit states without 
stabilizer measurements and to protect them both against lo- 
cal sources of decoherence and thermal fluctuations for times 
much longer than the longest time scale of the free qubit sys- 
tem, even when pulse imperfections are taken into account. 
Furthermore this method is fairly versatile, as it can be gen- 
eralized to qubit systems with Ising, XY , and Heisenberg in- 
teractions as well as to arbitrary (2D stabilizer) codes based 
on local stabiUzers in a straightforward manner. We provided 
simple, accurate analytical estimates for the gate and codeword 
fidelities of our method as a function of the system parameters, 
the period length T of the pulse sequence, and the magnitude 
of randomly distributed pulse errors. This allows us to esti- 
mate the maximum tolerable size of pulse errors and to opti- 
mize tunable system parameters, such as T or the coupling of 
the qubits to the cavity 5, in order to maximize the hfetime of 
codeword states. 
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Appendix A: Operations to Generate tlie Planar Code 

We present a particular choice of operations Ri and R2 that 
can be used to dynamically generate (one-quarter of) the pla- 
nar code when appUed to Hq as part of the sequence shown in 
Fig. 2(b). They are a generalization of the operations given 
in Sec. Ill to a quadratic qubit lattice of arbitrary size, i.e., 
their restriction to an interior unit cell equals Eqs. (6) and 
(7). The additional operations compared to the 2 x 2 lattice 
are needed to generate proper (third-order) boundary terms, 
while using strictly linewise rotations and entangling opera- 
tions. Schematic illustrations of these operations are provided 
in Fig. 7 for a finite lattice of 25 qubits. 

The generalized -Ri consists of nine steps, eight of which 
are shown in Figs. 7(a)-(h). Note that compared to Fig. 3, 
the lattice is rotated by 45° . After even rows of qubit terms are 
transformed from 2; to x by a 7r/2 rotation about the y axis, dis- 
joint pairs of qubits are entangled along one of the two sets of 
parallel diagonals [downward in Fig. 7(b)] so that pairs on dif- 
ferent diagonals are aligned along the perpendicular direction. 
The following three steps [Fig. 7(b)-7(e)] interchange one of 
the z and y terms along the perpendicular direction for every 



interior cell, while leaving the qubits that belong to bound- 
ary terms unchanged. In that context, it is important to notice, 
that one of the three qubits contributing to a boundary term 
does not actually he on that boundary. The rotations that in- 
terchange y and z terms [Fig. 7(d)] are performed about the 
X axis. Thus, by rotating the terms on the upper and lower 
boundaries to [Fig. 7(c)] prior to the interchange operation, 
their relative orientation is protected. After the interchange, 
the boundary qubits are rotated back [Fig. 7(e)]. The diago- 
nals, along which qubits are interchanged in step 4 are chosen 
in the following way: (i) they are perpendicular to the entan- 
gling bonds between the second-order terms, (ii) only one pair 
of qubits is interchanged per unit cell, and (iii) no y qubit be- 
longing to a boundary term is affected. This choice is always 
possible. The remaining pairs of qubits are entangled in step 
6 [Fig. 7(f)] to form third- and fourth-order terms. Finally, the 
qubits of those interior cells that were interchanged in step 4, 
have to be changed back [subfigure (h)] in step 8, again after 
protecting the relative orientation of the boundaries by a proper 
rotation |Fig. 7(g)]. The last step 9 (not shown) depends on 
whether plaque tte or star operators are about to be generated. 
All qubits on the same row have the same orientation after step 
8. The axis, about which each row has to be rotated in step 9, 
is therefore uniquely determined. 

Operation R2 eUminates unwanted single- and two-particle 
terms in the average Hamiltonian. For the Ri given above, it 
has to flip the signs of all qubits (i) on odd rows and (ii) on the 
intersections of even rows with those diagonals, whose z and y 
orientations were interchanged in step 4 of This operation 
can be decomposed into three steps, as illustrated in Figs. 7(i)- 
(k). First, odd rows are rotated from zio x [Fig. 7(i)]. A sub- 
sequent TT rotation [double arrows in Fig. 7(j)] about the x axis 
of all qubits along the diagonals of step 4 of R\ then flips the 
sign of qubits (ii). Finally, the same operation as in step 1 is 
performed once more on the qubits on odd rows [Fig. 7(k)]. 

The remaining three-quarters of the code Hamiltonian are 
readily given by 7r/2 rotations of all those steps about the cen- 
ter of the code. Although the choice of i?i and R2 is not unique 
and might by modified depending on a particular experimen- 
tal realization, the complexity of the operations cannot be de- 
creased by more than one or two steps. A simplification can 
be achieved, e.g., by choosing entangling operations and rota- 
tions that conform to the fourfold symmetry of the quadratic 
lattice. Since the resulting i?i and R2 are not as easily gener- 
alized to lattices of arbitrary size and not as easily appUcable 
in a linewise fashion, we do not consider them here. 



Appendix B: Second-Order Corrections to Average Hamiltonian 

The exact average Hamiltonian ilav = ^a"^ is given 

by an infinite expansion in orders of TA and can, in princi- 
ple, be evaluated to arbitrary orders using the Magnus expan- 
sion. For T <C this expansion converges quickly and only 
the lowest orders are relevant to describe the dynamics of the 
pulsed system. Furthermore, all odd-order contributions can 
be suppressed by symmetrizing the generating pulse sequence 
in time with respect to its center [32]. A symmetric sequence 
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FIG. 7. (Color online) (a)-(h) Schematic illustration of the first eight of nine steps of Ri to generate a quarter of the planar code Hamiltonian 
for a lattice of 25 qubits. When restricted to an interior unit cell (one is marked by a dashed rectangle) these steps reduce to those shown in Fig. 
3. (i)-(k) Decomposition of R2 into a sequence of linewise rotations. Single and double arrows represent rotations by n/2 and tt, respectively. 



for the planar code is obtained by reversing the order in which 
the different quarters are generated after every period T, as il- 
lustrated in Fig. 8. The corresponding Hamiltonians Qi with 
i — 1, ... ,4 are given by 



Q,^[l + S/A{b + b^)]Q,+cjob^b. 



(Bl) 



Taken as a whole, the resulting sequence consists of 32 op- 
erations, which give rise to 24 toggling frame Hamiltonians. 
With respect to time T, these Hamiltonians are mirror sym- 
metric (Hi — H25-i for 13 < z < 24), while the se- 
quence of 32 operations Rj is antisymmetric (Rj — R. 



33-j 



for 17 < j < 32). With this, the leading-order deviation from 
Eq. (9) is essentially given by H^y ^ (x (TA)^. However, as 
the sub-sequences to generate the (average Hamiltonians) Qi 
are symmetric as well, the corresponding second-order devia- 
tions scale as (T A /4) ^ and can be neglected for the calculation 
(2) 

of i?av ■ Hence, based on a sequence with effective toggling 
frame Hamiltonians Qi, we obtain 

H^? = ^ >^ (l-4i/2)[Qi,[Qfc,Q,]], (B2) 



384 ^ 

i<j<fc<;<8 



where S^i is the Kronecker delta and Qj>4 = Qg-j- Note that 
since the corresponding spin Hamiltonians Qj commute, all 



higher-order contributions to _ffav vanish identically for (5 = 
in which case Eq. (9) gives the exact average Hamiltonian. In 
the case of a finite coupling to the cavity, 
Eq. (10) using 



(2) 



evaluates to 



[Qi, [Qfc, Qj]] = -SiOo[26Qi+ujo{b + b^)]iQk -Qj). 



(B3) 



To derive an approximation for the decrease in the gate fi- 
delity F'^'^^ [see Eq. (11)] that is caused by the leading-order 



(0) 



correction -ffav , we transform the propagator exp[— 
^^av"*)] into the interaction picture with and h'^'^ playing 
the role of the noninteracting system and the interaction, re- 
spectively. The transformation of the latter can be evaluated 
analytically yielding 



H2\t) ^ A + cos{ujot)B + sm{LOQt)C 



(B4) 
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time 




-2T- 



<T/4> 



where R = R{0) is the perfect pulse operator. To estimate the 
statistical effect of these random errors on the fidelity of the 
whole planar code sequence, it is sufficient to consider the case 
of two inverse operations with independent errors flanking a 
period of free propagation: 



=R^{5e2)RiSei)R\5ei)Uo{t)R{S9i) 

{6e2)R{56r ) exp {-itR} {Sei)HoR{Sei )] 



(C2) 



FIG. 8. Symmetric pulse sequence for the generation of the planar 
code. The order in which the quarters Qi are generated is reversed 
after every period T. With respect to T, the 24 toggling frame Hamil- 
tonians Hi = TZlHoTZi are symmetric, while the 32 operations Rj 
are antisymmetric. 



The third term on the right-hand- side can be interpreted as a 
time evolution operator with effective Hamlltonian Hi + 6 Hi, 
where Hi = W {0)HoR{Q) is the perfect-pulse toggling 
frame Hamlltonian and 



with Hermitian operators 



A = 



B 



C 




3 



(B5) 



8A 



QV(6^ - b). 



To determine the approximate fidelity for times t ^ A~^, we 
replace the time -dependent second-order contribution by its 
average over the cavity's oscillation period Tc = 27r/a;o. It is 



given by m^' := T^^ ' m^'{t) dt = A. To second-order 
in A, the fidelity F^^^ can then be written as 



F(2)«|TrQB(l-ia-|A2)|/2^, 



(B6) 



where Tiqb denotes the trace over the qubit subsystem. Using 
TrQB(QfeQ() = 5kiL{L — l)/2 (the number of anyon opera- 
tors per code quarter for k = l),a straightforward calculation 
yields Eq. (12). 



Appendix C: Gate Fidelity in Presence of Pulse Errors 

To second-order in 66, an imperfect pulse R{S6) 
exp(i[^o + 56]S/2) can be approximated by 



RiSO)^[{l-'l)l+^'ls 



R, 



(CI) 



5Hi = --^{Hi-SHiS) + i^- [S, Hi]. 



(C3) 



In the case of Hq = —Aai^\ S — <y\l'' , and 9q = 7r/2, for 
example, application of R{59) leads to an average effective 



(1) 



Hamiltonian -A[(l - 56'^ /2)at' + dOa^^^]. The remaining 
terms from the right-hand- side of Eq. (C3) gives 



R\5e2)R{5ei) 

2^ 502 



i 



Sef + 2661662 + 66^ 



^ Ml- 662^ 



(C4) 



For the estimate of how much the statistical expectation value 
of the gate fidelity as a function of S6 deviates from the perfect- 
pulse value, we only need to consider terms that are propor- 
tional to the identity operator (as the trace of all PauU oper- 
ators vanishes) and 56^ (as for independent, unbiased errors 
{S61S62) = 0). The only such terms are contained in the 
first summand on the right-hand-side of Eq. (C4) and yield 
{F) « 1 — cTg /4 for a single pair of mutually inverse im- 
perfect operations. Consequently, for a nested sequence of n 
such operations, the deviation amounts to n times that value. 
Hence, by counting the number n{L) of mutually inverse oper- 
ations in the generating sequence for a planar code of length L, 
we obtain (F) w 1 — tn{L)(Tg/{4T) for stroboscopic times 
t = mT with integer m. Based on the operations shown in 
Appendix A, the number n{L) for the whole generating se- 
quence evaluates to n{L) = 36L^ — 20L + 14 and we obtain 
(F) « 1 — Cent — Cav^^ with Cgrr givcu by Eq. (13). 
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